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Abstract 

This paper demonstrates the progressive failure 
analysis capability in NASA Langley’s COMET-AR 
finite element analysis code on a large-scale built-up 
composite structure. A large-scale five stringer 
composite panel with a 7-in. long discrete source 
damage was analyzed from initial loading to final 
failure including the geometric and material 
nonlinearities. Predictions using different mesh sizes, 
different saw cut modeling approaches, and different 
failure criteria were performed and assessed. All 
failure predictions have a reasonably good correlation 
with the test result. 

Introduction 

For composite material to be successfully used in 
primary aircraft structures such as wings and 
fuselages, one of the design requirements is a two-bay 
crack damage tolerance capability. The two-bay 
crack, representing a discrete source damage caused 
by fragments from a failed engine, spans from one 
bay to the other bay and severs one stiffener. To meet 
the damage tolerance requirement, this two-bay 
damage cannot propagate through the adjacent 
stiffeners when the load level is below seventy 
percent of the design limit load. 1 * 

Analyses and tests have been conducted to 
investigate the damage progression and residual 
strength of stiffened composite panels with this type 
of damage. 1 ' However, robust analysis tools which 
can predict the material and geometrically nonlinear 
behavior of composite structures are still lacking. 
Good correlation between analytical results and 
experimental data for a large-scale stiffened 


composite structure with discrete source damage has 
yet to be achieved. 

In the past few years Moas 3 , Pifko 4 , Sleight 5 , and 
Sleight et al 6 implemented a robust damage 
progression analysis capability into the NASA 
Langley’s research finite element code, COMET- 
AR 7 . This code has an excellent nonlinear analysis 
capability and can perform adaptive mesh refinement. 
The damage progression analysis procedure, widely 
used by other researchers 7 14 , evaluates failure ply-by- 
ply and discounts the material properties for the failed 
ply. Other damage progression analysis approaches, 
such as the one based on continuum damage 
mechanics 1516 , can also be implemented. Using the 
COMET-AR’ s nonlinear solution procedure, the 
analysis can include both the material nonlinearity 
due to the degradation of material properties and the 
geometric nonlinearity due to the large deformations. 

Progressive damage analysis is computationally 
intensive. Failures are evaluated at each through-the- 
thickness integration point of every element, and 
material properties for each failed integration point 
need to be updated. Furthermore, to account for the 
geometrical nonlinearity and material degradation, the 
Newton-Raphson method 17 is used, and failure 
evaluations are performed at each converged load 
step to update the tangent stiffness matrix. Because of 
the complexity of the progressive failure analysis 
procedure for large-scale structures, it has rarely been 
successfully demonstrated and published. 

The objective of this study is to demonstrate that 
this newly implemented damage progression analysis 
procedure can be successfully used for analyzing a 
large-scale composite structure with discrete source 
damage. Analytical predictions of a stiffened panel 
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are compared with the panel’s test data to 
demonstrate that reasonable test and analysis 
correlation could be obtained. This panel was 
designed and fabricated by the McDonnell-Douglas 
Aerospace (MDA), now The Boeing Company, under 
NASA Langley Research Center’s “Innovative 
Composite Aircraft Primary Structure Program,” and 
tested at NASA Langley Research Center 1 . Another 
problem analyzed demonstrates that crack turning at 
the stiffener flanges could also be predicted. The 
crack turning at the flanges is a desired failure 
mechanism for an aircraft structure with a two-bay 
crack since the damage propagation can be confined 
within the two already damaged bays. 

A few technical issues, that are often encountered 
in the damage progression analysis of a large-scale 
stiffened composite structure with a discrete source 
damage, are also discussed. These technical issues are 
the effect of mesh size on the damage progression 
prediction, the effect of modeling the crack tips on the 
damage progression prediction, and the effect of 
using different failure criteria with their 
corresponding material degradation models on the 
damage progression prediction. 

By studying these technical issues, the usefulness 
of this progressive damage analysis method and its 
limitations can be better understood and defined. 

Failure Criteria 


Three failure criteria were used in this study 
including the Maximum Stress 18 , the Hashin 
criterion 19 , and the Hashin-Rotem criterion '. The 
Maximum Stress criterion is a limit criterion in which 
the stress interaction effect is not included for failure 
prediction. The Hashin and the Hashin-Rotem criteria 
maintain some stress interactions. These failure 
criteria can be expressed as follows; 
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Maximum Stress criterion 


Fiber failure: 
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Shear matrix cracking: 


The principal stresses and strength properties of a 

lamina used in the three failure criteria are defined as 

o n , o 2 t , and x 12 : stress in the fiber direction, stress 
in the transverse-to-the fiber 
direction, and the in-plane shear 
stress, respectively. 

X and X’ : tensile and compressive strengths in fiber 
direction. 

Y and Y' : tensile and compressive strengths in 
transverse direction. 
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S A and S T : in-plane shear strength and transverse 
shear strength, respectively. 

All three-failure criteria have the capability of 
identifying failure modes of a lamina. In the 
progressive damage analysis, the identification of 
failure modes in a lamina is critical for the 
continuation of progressive analysis. Failure modes 
will be used to determine which material properties of 
the failed lamina need to be degraded for the next step 
of the analysis. Many other failure criteria 21 22 can also 
used for the damage progression analysis. 


such that the stiffeners and flanges were not directly 
loaded. Dimensions of the panel are shown in Figure 
2. The length and the width of the panel are 80 and 40 
inches, respectively. The height of the stiffeners is 
2.30 in. and the width of the flanges is 3.2 in. The 
stringer spacing is 8.0 in. This panel has five-stringers 
and a 7-in. saw-cut to simulate a discrete source 
damage. The saw-cut severs the center stringer and 
3.5 in. of skin on either side of the stringer centerline. 
In the analysis, moduli and strengths for a typical 
AS4/3501-6 unidirectional lamina were used. They 
are: 


Material Degradation Models 

The material degradation model used in this study 
for each aforementioned failure criterion is listed in 
Table 1. For example, when fiber failure is detected, 
the moduli E n and G 12 and the Poisson’s ratio v 12 
are degraded according to the Maximum Stress 
Criterion. Similarly, when matrix failure or shear 
failure are detected, the moduli E 22 and G 12 and the 
Poisson’s ratio v 12 are degraded. 

The material properties can be degraded slowly or 
can be instantaneously reduced to zero if failure 
occurs. In this study, the ^1 1 ' ^22 > and G p were 
instantaneously reduced to 1,000 psi and v 12 was 
instantaneously reduced to 0.01. Retaining small 
values of material properties was found to help avoid 
numerical difficulties in the nonlinear analysis 
solution. 


E x = 15.04 msi 
Et = 1 .60 msi 
G 12 = 0.8 msi 
v 12 =0.34 

Ply Thickness = 0.0064 in. 


X = 182.0 ksi 
X' = -138.97 ksi 
Y =9.94 ksi 
Y' = -31.04 ksi 
S A = S T =17.52 ksi 


The lay up sequences are [45/-45/0 2 /90/0 2 /-45/45] 6 
for the skin, [45/-45/0 2 /90/0,/-45/45] s for the stiffener, 
and [45/-45/0 2 /90/0 2 /-45/45] 4 for the stiffener flange. 
The panel was fabricated by a stitched/RFI (resin film 
infusion) process. Because the stitching can prevent 
the delamination failure mode, the present analysis 
procedure, which did not evaluate the delamination 
failure mode, was considered to be appropriate for 
predicting the damage progression in the stitched 
panel. 


Nonlinear Analysis 


Failure in an element is quantified as the percent 
of ply-integration point failure in that finite element. 
The COMET-AR element used for this analysis is 
ES5/E410, which is a 4-node quadrilateral element. 
Each element has five Gauss quadrature points. 
Through-the-thickness integrations of the element use 
three integration points for each ply in the laminate. 
The total ply-integration points in an element are 
5x3xNumber-of-Plies. In this paper, a damaged 
element means that more than 60% of its ply- 
integration points have failed, resulting in a very 
small load carrying capability of this element. 

Panel Description 

The MDA panel in testing is shown in Figure 1. 
The applied load was introduced to the panel through 
two metallic test fixtures, one at the top and one at the 
bottom. The test fixtures were bolted to the skin bays 


The nonlinear behavior of a composite structure 
can be caused by geometrical nonlinearity, material 
nonlinearity, or a combination of both. Geometrical 
nonlinearity is associated with a large deformation 23 . 
Material nonlinearity can occur when the stress-strain 
relations become nonlinear due to ply failure. In the 
progressive failure analysis of a composite structure 
with a large deformation, both the geometrical 
nonlinearity and material nonlinearity need to be 
included. 

The assembled nonlinear finite element system 
equation is 

[x,({d})Xd}=H (11) 

where {/)} is the displacement vector, {p} is the 
applied load vector, and [x, ({/)})] is the tangent 
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stiffness matrix which depends on the nodal 
displacements. 

The Newton-Raphson method l : was used to 
solve this nonlinear system equation (Equation 11). 
The damage evaluation and the recomputation of 
\K t \ were performed at every converged load step. 
Equilibrium between the internal force vector, left 
side of the system equation which was computed by 
using the new [/if,], and the external force vector, 
right side of the system equation, was reestablished 
by repeating the nonlinear analysis and failure 
evaluation at the current load step. After the 
equilibrium was established, the nonlinear analysis 
advanced to the next load step. 

Finite Element Models 

Three finite element models were created. They 
are a coarse mesh crack model, a fine mesh crack 
model, and a fine mesh slit model. The three models 
are shown in Figures 3 to 5, respectively. These 
models represent the fourth quadrant of the panel 
shown in Figure 2. In the first two models, the width 
of the saw-cut was neglected and the discrete source 
damage was modeled as a sharp crack. In the third 
model, the width of the saw-cut was modeled as a 
0.2-in. wide slit with round tips. The radius of the 
round tip is 0.1 in. Boundary conditions used in the 
analyses are shown in Figures 3 to 5. Symmetry 
boundary conditions were applied at the uncracked 
portion of the top edge and the right edge, free 
boundary conditions were applied at the left edge, and 
a uniform longitudinal displacement load was applied 
at the bottom edge of these finite element models. 
Note that loads were applied on the skin bays only 
and not applied on the blades and the flanges to 
simulate the actual test loading conditions. 

In the region surrounding the crack, the mesh size 
of the fine mesh crack model is half the mesh size of 
the coarse mesh crack model. The width of the coarse 
mesh is 0.2 in. while the width of the fine mesh is 0.1 
in. For the fine mesh slit model with a round tip, the 
mesh size at the tip is much less than the 0. 1 in. while 
the mesh size away from the crack tip is the same as 
the fine mesh crack model. The coarse mesh crack 
model contains 1,302 nodes and 1,238 quad element, 
the fine mesh crack model contains 3,248 nodes and 
3,163 quad elements, and the fine mesh slit model 
contains 3,348 nodes and 3,254 quad elements. 


Results from the coarse mesh crack model and the 
fine mesh crack model were used to study the mesh 
size effect on the damage progression. The issue of 
modeling the saw cut as a sharp crack or as a slit with 
round tips was studied by using the two fine mesh 
models. The issue of how different failure criteria and 
their corresponding material property degradation 
models could affect the damage progression results 
was studied by using the results from all three 
models. The Maximum Stress, Hashin, and Hashin- 
Rotem criteria and their corresponding material 
property degradation models were used in the 
analyses. To demonstrate the potential of using this 
damage progression analysis capability, results 
predicted with these three models were compared 
with the MDA panel test results to illustrate that a 
reasonably good analysis and test correlation can be 
obtained. To demonstrate that crack turning at the 
flange can be predicted by using this progressive 
failure analysis capability, the fine mesh crack model 
was used for the crack turning study. 

Damage Progression Analysis Results 

During the nonlinear analysis, a longitudinal 
displacement was applied on the bottom edge of the 
panel’s skin bays in a small increment of 0.002 in. 
The applied load level, that is the reaction load for 
each load step, was computed. Applied load levels for 
the advancing damage to reach the flange predicted 
by three different finite element models and three 
different failure criteria are reported in Table 2. A 
typical damage growth plot is shown in Figure 6. 
Note that a damaged element means that more than 
60% of its ply-integration points failed. The damage 
growth can also be visualized from a deformed plot as 
shown in Figure 7. The distortion and the large 
deformation of the damaged elements ahead of the 
crack can be observed. The deformed shape shown in 
Figure 7 was obtained by mirroring the deformed 
shape of a quarter model. Results for the Figures 6 
and 7 were obtained by analyzing the fine mesh 
model with the Hashin-Rotem failure criterion. 
Degradation of material properties of the damaged 
elements may result in an ill-conditioned tangent 
stiffness matrix, and large deformations of these 
damaged elements may also cause divergence 
problems in the nonlinear analysis. An approach used 
by Chamis et al 1 ’ to adaptively remove the severely 
damaged elements from the model may prevent the 
divergence problem. This approach will be 
incorporated into COMET-AR code in the future. 
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Mesh size effect 

Table 2 shows that the applied load levels 
predicted by the coarse mesh crack model are higher 
than those predicted by the fine mesh crack model. 
Since the high stress field exists at the crack tip, the 
finer mesh predicted higher element stresses, which 
resulted in an earlier damage growth initiation in the 
fine mesh crack model. 

As illustrated in Figure 8, a slower damage 
growth for the coarse crack model was found. At a 
load level about 485 kips, the damage growth for the 
fine mesh was 1.1 in. while the damage growth for 
the coarse mesh was 0.8 in. Note the damage growth 
length is defined as the length of a damage zone 
measured from the original crack tip along the crack 
line. For the damage propagating to the edge of the 
flange of the adjacent stiffener as shown in Figure 9, 
the applied loads for the coarse mesh crack model and 
the fine mesh crack model were 588 and 564 kips, 
respectively. A 4.3% higher load was needed for the 
coarse mesh crack model. Results in Figures 8 and 9 
were obtained by using the Hashin-Rotem failure 
criterion. 

Crack tip modeling 

Results in the second row and the third row of 
Table 2 show how the crack tip modeled can affect 
the damage progression. The damage progression 
results obtained from this fine mesh slit model were 
compared with the results obtained from the fine 
mesh crack model. The applied load levels predicted 
by the fine mesh crack model were higher than the 
fine mesh slit model with a round tip. The difference 
may be caused by a much earlier damage initiation 
for the fine mesh slit model due to a very fine mesh 
used to model the round tip. 

As illustrated in Figure 10, at a load level about 
360 kips; the damage growth length for the fine mesh 
slit model was 0.6 in. while for the fine mesh crack 
model it was 0.2 in. For the damage progression to 
the edge of the adjacent stiffener as shown in Figure 
11, the applied loads for the fine mesh slit model and 
the fine mesh crack model the were 522 and 564 kips, 
respectively. An 8.0% higher load was needed for the 
fine mesh crack model. Figures 10 and 11 were 
obtained by using the Hashin-Rotem failure criterion. 


Effect of failure criterion used 

Results in Table 2 show that the damage growth 
depends on the failure criterion and its corresponding 
degradation model used. Comparing the damage 
progression results predicted by the Maximum Stress 
criterion, the Hashin criterion and the Hashin-Rotem 
criterion, the highest load level for every mesh model 
used was predicted by the Maximum Stress failure 
criterion. The lowest load level for every mesh model 
used was predicted by the Hashin failure criterion. As 
mentioned before, the Maximum Stress failure 
criterion is a limit criterion in which the failure of a 
ply is predicted by comparing the principal stresses in 
that ply with its principal strengths, separately. In the 
Hashin criterion and the Hashin-Rotem criterion, 
some stress interactions are maintained. The shear 
stress can contribute to the fiber and matrix failure for 
the Hashin failure criterion while the shear stress can 
only contribute to matrix failure for the Hashin- 
Rotem criterion. 

As illustrated in Figure 12, the load levels 
predicted by using the fine mesh slit model with the 
Hashin and the Maximum Stress criteria were 487 
kips and 547 kips, respectively. The difference is 
about 12.3%. For the same model, the load level 
predicted by the Hashin-Rotem was 522 kips as 
shown in Figure 11, which is 4.6% below the applied 
load level predicted by the Maximum Stress failure 
criterion. 

Comparison with test result 

The test result documented in Reference 1 shows 
that at a load level of 521 kips (96% of design limit 
load), the damage progression reached to the flange 
of the adjacent stiffeners. Note that this five-stringer 
panel meets the discrete source damage requirement, 
which requires the damage progression confined 
within the two damage bays for the applied load 
below 380 kips (70% of the design limit load). 

As discussed in the previous three subsections, the 
mesh size, the different modeling approach of the 
saw-cut, and the selection of the failure criterion can 
have some effects on the predicted damage growth 
results. In general, the applied loads predicted by all 
the models using the Maximum Stress, the Hashin 
and the Hashin-Rotem criteria are considered within 
an acceptable range from the test result. Higher loads 
were predicted when the coarse mesh crack model or 
the Maximum Stress criterion was used. The worst 
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prediction is about 16.3% above the test result when 
the coarse mesh model and the Maximum Stress 
criterion were used. The best prediction is about 0.2% 
above the test result when the fine mesh slit model 
and the Hashin-Rotem criterion were used. 

Potential for crack turning prediction 

For a stiffened composite panel to meet the two- 
bay damage tolerance requirement, the crack growth 
needs to be contained within the two damaged bays 
and can not propagate through the adjacent flanges. It 
is desired that the crack turns when it reaches the 
adjacent stiffener flanges rather than penetrating the 
flanges. It is expected that if a high strength material 
system is used in the stiffener flanges, the crack can 
be stopped and forced to propagate parallel to the 
stiffener. In this study, an assumed high strength 
material system was used for the flange. Indeed, the 
crack changes direction at the edge of the flange 
region as shown in Figure 13. In this study the 
Hashin-Rotem failure criterion was used. 

Concluding Remarks 

Analytical results documented in this paper 
demonstrate that the recently implemented 
progressive failure analysis capability in COMET-AR 
can be a useful tool for the damage tolerance analysis 
of large-scale built-up composite structures. A large- 
scale five stringer composite panel with a 7-in. long 
discrete source damage was analyzed. Geometrical 
and material nonlinearities were considered in the 
analysis. Predictions using different mesh sizes, 
different saw cut modeling approaches, and different 
failure criteria were performed. All the predictions 
have a reasonably good correlation with the test 
result. 

A few technical issues, which can affect the 
damage progression predictions, were also discussed. 
Analysts need to be aware of these issues and 
exercise the progressive damage growth analyses with 
care. 

The first technical issue is the effect of the mesh 
size. A coarse mesh crack model and a fine mesh 
crack model were used to study the mesh size effect. 
It is found that for the damage to propagate to the 
same distance, a higher applied load was needed for 
coarse mesh model. 


The second technical issue is the effect of crack 
tip model. The discrete source damage can be 
modeled as a sharp crack or as a slit with round tips. 
A fine mesh crack model and a fine mesh slit model 
were analyzed. Results show that the damage 
initiation is much earlier in the slit model resulting in 
a longer damage growth when both models were 
subjected to the same load level. This indicates that 
the configuration of the saw-cut needs to be carefully 
modeled. 

The third technical issue is the effect of different 
failure criteria. Three different failure criteria were 
used to study how the damage progression can be 
affected by the selection of failure criterion. For the 
same amount of damage growth, the applied load 
level predicted by the Maximum Stress criterion is 
higher than that predicted by the Hashin and the 
Hashin-Rotem criteria. The Hashin-Rotem criterion 
provided the closest match with the experimental 
results. 

Further improvement of the accuracy of 
predictions is needed before analysis simulations can 
be used to replace costly structural tests. This may be 
achieved by the advancement of micromechanics to 
better understand the failure mechanisms, more 
advanced failure analysis procedures, and 
incorporating reliability analysis in the failure 
prediction. Nevertheless, this progressive failure 
analysis capability should be able to provide pretest 
predictions for damage tolerance tests of aircraft 
structures. These predictions are vital for the selection 
of a test machine, fixture design, and the set-up of test 
data acquisition system. 
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Table 1. Material Degradation Models 


Failure Criterion 

Fiber failure 

Matrix failure 

Shear failure 

Maximum Stress 

E n ' ^12 ’ V 12 

E 22 > ^12 ’ V 12 

E 22 ’ G\ 2 > V 12 

Hashin 

E n > G u , v 12 

E 22 , G n 7 v 12 

Not applicable 

Hashin-Rotem 

E u ’ G\2 >V[2 

E 2 2 > G l2 , v 12 

Not applicable 


Table 2. Applied Loads for Damage Progression to Reach the Adjacent Flange 


FEM Models 

Maximum Stress 
Criterion 

Hashin Criterion 

Hashin-Rotem Criterion 

Coarse mesh crack 
model 

606 kips*(16.3%)** 

571 kips (9.6%) 

588 kips(12.9%) 

Fine mesh crack model 

591 kips (13.4%) 

553 kips (6.1%) 

564 kips(8.3%) 

Fine mesh slit model 

547 kips (5.0%) 

487 kips (-6.5%) 

522 kips(0.2%) 


* Thousand pounds. 

** % difference relative to the test data of 521 kips. 



Figure 1. A five-stringer composite panel with a discrete source damage in testing. 
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Figure 2. Geometry of the five-stringer composite panel. 


V = 0.0 



(BC’s on the bottom edge were applied on the skin bays only) 
Figure 3. The coarse mesh crack model. 
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V = 0.0 



Figure 4. The fine mesh crack model. 
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(BC’s on the bottom edge were applied on the skin bays only) 


Figure 5. The fine mesh slit model with a round tip. 
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Figure 6. Damage reaching the edge of the adjacent flange at a load level of 564 kips. 



Figure 7. Deformed shape and local damage zone. 
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Figure 8. Damage growth advancing faster for the fine mesh model. 
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Edge of Flange 




Figure 9. Applied load level predicted by the coarse mesh higher than that predicted by the fine mesh model. 
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Fine mesh slit model, 363 kips Fine mesh crack model, 369 kips 


Figure 10. Damage growth advancing faster in the fine slit model at the initiation stage. 
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Fine mesh slit model, 522 kips 


MMiflSt 

Fine mesh crack model, 564 kips 


Figure 1 1 . Applied load level predicted by the fine slit model lower than the fine crack model. 
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Edge of Flange 



Hashin criterion, 487 kips 


Maximum stress criterion, 547 kips 


Figure 12. Applied load level predicted by the Hashin criterion lower than that by the Maximum Stress criterion. 
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Figure 13. Crack turning due to using high strength material in the flange. 
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